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The effects of graphite surface geometrical deformation on the dynamics of conducting 
electrons are investigated theoretically. The analysis is performed within the framework of a 
deformation-induced gauge field and corresponding deformation-induced magnetic field. It 
is shown that the latter gives a local energy gap along the axis of a deformed nanotube. We 
compare our energy gap results with experimental data on energy gaps in nanotubes and 
peapods. We also discuss the mixing of two Fermi points and construct a general model of 
low energy dynamics, including a short-range deformation of the graphite sheet. This model 
is equivalent to the Weyl equation in U(l) Abelian and SU(2) non-Abelian deformation- 
induced gauge fields. 

§1. Introduction 

The geometric structure of materials and their electronic and magnetic prop- 
erties are closely related. It is widely recognized that carbon nanotubes 1 ) provide 
us with a great opportunity to study this relationship. For instance, a single-wall 
carbon nanotube (SWNT) exhibits either metallic or semiconducting behavior, de- 
pending on the lattice structure around the tube axis. 2 ) This structure dependent 
electric property can be understood by noting that the energy band of an SWNT is 
determined by the wave vector quantized by the periodic boundary condition around 
the tubule axis and the energy dispersion relation of graphene. 3 ) 

The energy dispersion relation of a plane of graphene is unique, because of the 
linear dispersion relation at the two distinct Fermi points (the K and K' points). 
Moreover, since the unit cell of a honeycomb lattice structure consists of two sublat- 
tices, the linear energy dispersion relation leads to dynamics described by the Weyl 
equation. 4 ) 

However, the dynamics of the conducting electrons of realistic graphene-based 
materials are different from those of flat graphene, because in the former case, 
there are shape fluctuations that result in the modification of the overlap integral 
of nearest-neighbor 7r-orbitals. Even when a graphene sheet is folded into a cylinder 
to form a carbon nanotube, the mean curvature of the graphite surface changes the 
dynamics of the conducting electrons. 5 ) In this perturbation due to atomic 

deformation appears in the dynamics through the gauge coupling. 6 ) A conducting 
electron moves in the curved graphite sheet in which the effect of curvature can be 
expressed in terms of a gauge field. The effective gauge field shifts the Fermi point 
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in the Brillouin zone due to the Aharonov-Bohm effect and therefore changes the 
energy band structure. The effect of the finite curvature on the low-energy dynamics 
was elucidated in this way by Kane and Mele. 6 ^ Their theoretical prediction for the 
energy gap induced by the curvature effect was confirmed by the scanning tunneling 
spectroscopy (STS) experiment of Ouyang et al. 7 ) and was applied to different chiral 
structures in several subsequent works. 8 ) 

In this paper, we generalize the idea of the effective gauge field to a local defor- 
mation of the graphite surface and show that the deformation-induced gauge field 
for a locally deformed lattice provides a local energy gap. This gauge field is useful 
for describing finite-scaled nano-meter materials. 

This paper is organized as follows. In §§|2]and|2J we illustrate how the electron 
dynamics depend on the geometry of a deformed surface. In particular, we consider 
how a local deformation affects the local electrical properties of conducting electrons 
in SWNTs in terms of the deformation-induced gauge field. We point out here that 
the deformation-induced gauge field and the usual electro-magnetic gauge field are 
different, though the analogy is useful in many cases. Furthermore, the deformation- 
induced magnetic field concerns not the spin of an electron but its pseudo-spin, 9 ) 
as defined by the two sublattices of graphite and SWNTs. In § |IJ we clarify the 
relationship between the atomic deformation and a local energy gap structure. We 
compare our theoretical result with the local energy gap observed by STS in nan- 
otubes 7 ) and peapods. 10 ^ In § |3 we consider a short-range deformation that results 
in a mixing of the wave functions at two Fermi points, and we construct an effective 
model for the low-energy dynamics. 

§2. Definition of the Gauge field for deformation 

We consider the quantum behavior of conducting electrons on a graphite surface 
using the nearest-neighbor tight-binding Hamiltonian, 

Wnear = ^ ^ V a (ri)a\ +a ai + h.c. (1) 

Here, A (in the summation index) denotes an A-sublattice, and a\ are the canon- 
ical annihilation-creation operators of the electrons at site i, which satisfy the anti- 
commutation relation {aj,aj} = 5{j, and a site labeled i + a (with a = 1,2,3) is 
a nearest neighbor of site i. We include deformation of the graphite surface in the 
form of the position-dependent hopping integral V a (r^) in the Hamiltonian. 

We decompose the hopping integral into two components as V a {ri) = V n +5V a (ri) 
and define 7i ne ar = 7~Lo + Hdefomu where 

Ho = X)E^°!+« o *+ h - c -- ( 2 ) 

a ieA 

^deform = ^ ^ 5V a (ri)a\ +a ai + h.c. (3) 

Hereafter, we refer to the deformed Hamiltonian as "Hdeform- The total Hamiltonian 
adopted for this study is TCneur = Wo + Wdeformj which is not easy to solve for general 
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8Va(?i)- However, as long as we consider behavior near the Fermi level, we can obtain 
several important physical consequences; for example, the model Hamiltonian can 
exhibit a local energy gap 10 ) for a specific deformation 5V a (ri). 

In this paper, we mainly consider the case in which the deformation is sufficiently 
delocalized in comparison with the diameter of an SWNT, so that there is no inter- 
action between two Fermi points. In this case, we have doubly degenerate electronic 
states near the K and K' points. We follow the effective-mass description 11 ) for 
low-energy conducting electrons around each Fermi point. Then, the Hamiltonian 
around the K-point, Hk, is given by 

H K = v F a -(p-A), (4) 

where a ■ p = u\p\ + 02P2 and a ■ A = o\A\ + a 2 A 2 . Here, A is a vector field defined 
on the surface, up is the Fermi velocity and a% (i = 1,2,3) are the Pauli matrices, 
given by 

ai = (J J) ' CT2 = (J 7) ' °" 3 = (J -l) • (5) 

We define coordinate axes around and along the nanotube axis as x\ and x%. The 
first term in Eq. Q, vpo -p, describes the energy dispersion of Tio, 4 ^ and the second 
term, —vpo ■ A, denotes the energy dispersion of Wdeform- The vector field vfA(t) 
is a linear function of SV a (r). 

We derive Eq. Q and obtain an explicit form of A(r) in the next section, but here 
we note two points. First, the deformation can be included in the Hamiltonian as a 
gauge field, because of the gauge coupling p — ► p — A. We call A the deformation- 
induced gauge field and distinguish it from the electro-magnetic gauge field A era . 
Note that A does not break the time-reversal symmetry, because the sign in front of 
A is different for the K and K' points. This contrasts with the fact that A era does 
violate the time-reversal symmetry (see Eq. I22|) . The second point is that we can 
decompose any gauge field into the sum of constant, rotationless, and divergenceless 
components as 

Ai = A , + d^a + etjdjVi, (i, j = 1,2) (6) 

where is the antisymmetric tensor (ei2 = —£21 = 1 and en = 0), and and 3^ 
are regular scalar functions defined on the tube surface. 

Several important physical properties can be derived from through dimen- 
sional reduction. We integrate TLk over the circumferential coordinate (xi) to obtain 
a one-dimensional theory: 

j> ^J^K = vpipl - A\- 5m(x 2 ))a 1 + v F (p 2 - A° 2 )a 2 . (7) 

Here, |C^| is the circumference of the nanotube, p\ is determined by the chiral index, 
and we have 

/dx\ 
^r^d2%{x l ,x 2 ). (8) 
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The function H/ a in A\ disappears in the integration over the xi coordinate. We have 
ignored the \P a part of A2, because it can be expressed as the phase of the wave 
function. 

We identify the o\ term with the mass of a particle propagating along the axis 
direction (a^). This is because (1) the energy spectrum at position X2 is bounded by 
the square root of the sum of the squares of each coefficient of the Pauli matrices in 
Eq. as ±vpy/ (p® — A® — 5m{x2)) 2 + (p2 — A®) 2 , and (2) p2 scales as the inverse 
of the system length (1/|T|), and hence P2 — A® can be ignored in comparison with 
Pi — A\ — 5m(x2) for a long system (|T| 3> |Ch|). We therefore define the local energy 
gap as the minimum energy gap for long tubes, given by 

£gap(z 2 ) = 2v F \p\ — A\ — Sm(x 2 )\ = 2v F 

Because the circumferential component (Ai) of the deformation- induced gauge field 
is determined by the modulation part of the hopping integral, knowing 5V a (r) allows 
us to calculate the local energy gap using the above formula. The modulation part of 
the hopping integral can be calculated by choosing a specific scheme. We adopt the 
Slater-Koster scheme and apply the formula to several deformed nanotubes. In the 
next section, we give a detailed derivation and present a general formula applicable 
to an arbitrary chiral index. 

§3. Derivation of the local energy gap 

There are two main purposes of this section. One is to obtain an explicit form 
of A in terms of 6V a (ri), which is necessary to derive an explicit formula for the 
local energy gap, and the other is to state a general property of the low-energy 
Hamiltonian. 

For the above-stated purposes, we first rewrite Eq. (j2J) using the Bloch basis 
vectors as 

Ho = V 7r J2f( k )\^A}^B\+h.c, (10) 

k 

where we define f(k) = J2 a fa{k) an d f a {k) = e lk ' Ra . Here, R a (a = 1,2,3) are 

vectors pointing to the nearest-neighbor sites from an A site (Fig. We have 
defined the Bloch basis vectors with wave vector k as 

|**> = -^=E e * r4a » t l°>' t 11 ) 

V a iea 

where the subscript a (= A, B) denotes two sublattices, as illustrated in Fig. ^ and 
./V^ (= Nb) denotes the number of black (white) sites represented by solid (empty) 
circles. Here, T\ labels the vector pointing toward the site i. 
We expand f(k) in Eq. Q1U|) around the K point, kp, as 



Pi 



dx 

Jc, 



-Ai(xi,x 2 ) 



(9) 



f{k) = Y,fa(k F )i{k-k F )-R a + -- - , 

a 



(12) 
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Fig. 1. Structure of a honeycomb lattice with the two symmetry translation vectors cti = V3a cc e x 
and a2 = (\/3/2)a cc e x + (3/2)a cc e y , where e x and e y are unit vectors, and a cc is the nearest- 
neighbor bond length. The black (white) circles indicate the A (B) sublattices. The vectors R a 
(a = 1,2,3) point to the nearest-neighbor sites of an A site. They are given by _Ri = a cc e y , 
R2 = — (V3/2)a cc e x — (l/2)a cc e y and R3 = (V3/2)a cc e x — (l/2)a cc e y . 



where we use the condition for the Fermi point /(fc^) = 0. [We can set fi(kp) = 1, 
.M^f) = e~ 1 ^ and fzikp) = e +l ~%~ .] The lattice structure of SWNTs is specified by 
the chiral and translational vectors defined by Ch = na\ + ma^ and T = pa\ + qa2, 
where a\ and 02 are translational symmetry vectors in the planar honeycomb lattice 
(Fig. P). The corresponding wavevector k can be decomposed as k = ^i\k\ + \i2k21 
where [i\ and \ii are integer coefficients of the vectors k\ and k2, which satisfy 
Ch-k\ = 2-7T, Ch ■ &2 = 0, T ■ ki = and T ■ k2 = 2tt. With these definitions, the first 
term on the right-hand side of Eq. (jl()j) is approximated by 



V n e- l9 w\^)(^l 

where w = w\ — iu)2, Wi and W2 are given by 
3a, 



2\C h \ 

and the angle 9 satisfies 



(2717/1 - k F ■ C h ), w 2 



3a cc 
2\T\ 



(2-K(i 2 - k F • T), 



(13) 



(14) 



. rn\ .VS I C h | _ 

n-\ — 1 — m = —= — e 

2 J 2 V3a cc 



(15) 



To obtain Eq. (|13|) . it is useful to rewrite R a in terms of the chiral and translational 
vectors as 



Ri = 


3N C i 


R 2 = 


2 

3N C [ 


R 3 = 


2 - 
3N C l 



p+ V Ch 



m . 
n+-)l 



P - +q )c h+ (l +m] T 



(16) 
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where we introduce N c = mp—nq and use the relationship TxCh = (SyE/2)a^ c N c (e x x 
e y ). Finally, the leading term in the Hamiltonian matrix for TIq can be expressed in 
the basis of the two sublattices *(iZ%^b) as 

Next, we consider the momentum transfer matrix element of Wdeformj 

(* k A +Sk \n dc{orm \v k B ) = ^EE^(^)/aWe-^ r s (18) 

where we have used Eqs. (j3J) and (|lll) to obtain the right-hand side. Here, we restrict 
ourselves to considering the momentum transfer matrix element that does not mix 
the two Fermi points (i.e. we assume \8k\ <C |&f|). We then set k = kp + 8k' and 
obtain 

(* k A +Sk \H dc{oTm \# k B ) = ^ E (E<^H/a(M) e~ m - r *, (19) 

where we omit an O {8V a 8k') correction. The momenta k and k + 8k in the last 
equation should be recognized as being restricted to the wave vector around the K 
point. Up to leading order, this shows that the electrons near the K point are subject 
to the "local potential" given by ^2 a 8V a {ri)f a (kp), because its Fourier component 
appears in the right-hand side of Eq. (|19|) . It is important to note that there is 
a similarity between the deformation Hamiltonian and the potential V{x) (whose 
constant component should be subtracted) for a particle of mass m. Let us consider 

2 

a simple one-dimensional Hamiltonian, H = + V . When we consider the matrix 
element of the potential term between free particle states denoted |A:) with a definite 
momentum, we obtain (k + Sk\V\k) = J ^-V(x)e~ lSk ' x , where L is the system length. 
This should be compared with Eq. (|T§)) . and it shows that ^ a 8V a (ri)f a (kp) acts as 
a potential. Note also that Eq. ()17|) corresponds to the kinetic term of Eq. ()1U|) . 
Combining Eq. (|17|) with Eq. (|19|) . we obtain the effective-mass Hamiltonian 

^ f e- i9 w\ 



e w 
e~ ie 



i6 V F 



a ■ P - ( o e -i6 ) ° ■ A 



(20) 



where we have denoted ^ a 8V a (ri)f a (kF) as —vp(A x — iA y ) and used /i(&f) = 1, 
fiikp) = e~ l ~ and fz(kp) = e +% ~ . Here, we have defined the Fermi velocity as 
vf = 3V 7T a cc /2h, and we have a-p = (2h/3a cc )(aiWi+a2W2) and a- A = o\A x -\-<xiAy. 
Then, functions A x and A y are given by 



-(',.-,!,. - 8\\ - X -W 2 - ±8V 3 , —vpAy = ^ ( ( M - 8\}o. (21 ) 
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This expresses the explicit relationship between the hopping integral modulation 
and the deformation-induced gauge field. We refer to functions A x and A y , or their 
appropriate linear combination, as the deformation-induced gauge field. This is 
because the effects of the deformation can be included in the theory through the 
substitution p — > p — A. This is the same substitution as for an electro-magnetic 
gauge field. 

We simplify the notation in Eq. ()2()ll and state a general property of the effective 
Hamiltonian. We have two distinct modes, whose dynamics are approximated by the 
following two effective Hamiltonians in the presence of an external electro-magnetic 
gauge field A em : 

n K = v F a-(p-A-A era ), Hk' = vpo' ■ (p + A — A em ) . (22) 

Here, we have defined a = (o"i,<T2) and cr' = (— 0"i,C2)- The Hamiltonian Hk' can 
be obtained by expanding Eq. (JTU|) around the K' point, —kp, and repeating the 
calculation presented above. Our notation for the momentum is a ■ p = o\p\ + 02P2 
and a' ■ p = —o\p\ + 02P21 where 

H(2nm - k F ■ C h ) h{2irfx 2 — kp ■ T) 
Pi = j^j , P2 = jyj , (23) 

and our notation of the deformation-induced gauge field is a ■ A = a\A\ + 02A2 and 
a' ■ A = —u\A\ + (J2A2, where 

AA ( cos 9 sm9\ (A x \ 

A 2 J \-sm6 cos e)\Ay )' { > 

Note that p\ and P2 are, respectively, the momenta around and along the axis as 
measured from the K point. For the K' point, kp and [i\ should change sign, and, 
correspondingly, p\ for the K' point will have the sign opposite to that of the K 
point. Note also that the sign of A is different for the K and K' points, in contrast 
to the case for A em in Eq. (|22j) . Therefore, in the absence of an external electro- 
magnetic gauge field, time-reversal symmetry is preserved even in the presence of a 
deformation. 

By virtue of the decomposition given in Eq. @ with the gauge coupling structure 
of Eq. I)22|) . the following two points become clear. First, the \P a component does not 
change the energy spectra of the theory, because it can be eliminated by multiplying 
the wave function by a phase (i.e., taking tp — ► e~^tp). However, it should be 
mentioned that a nontrivial functional form of gives a nonvanishing divergence 
of the deformation-induced gauge field, represented by V • A = V 2 !^ a 7^ 0. This may 
correspond to the deformation potential and could cause a local (background) charge 
modulation. 12 ^ The effective one-dimensional quantum field theory 13 ) derived from 
Eq. (|22|) can be used to study this effect. Analysis using this theory shows that the 
'Fa component does not change the energy spectrum of the theory but gives only 
a charge density modulation for metallic energy band structure cases. Second, the 
function operates as a source of a deformation-induced magnetic field and can lead 
to an important physical effect. The curl part of the deformation-induced gauge field 
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defines a deformation-induced magnetic field as 

B± = etjdiAj = V 2 ^, (25) 

whose direction is perpendicular to the graphite surface, and a non-trivial magnetic 
field gives a local modulation of the energy gap, as defined in Eqs. (jSJ) and ©• 

In the previous section, in Eq. (J7|) we integrated over the circumferential coor- 
dinate (xi) of Eq. (|2*2*|) . In Eq. (JJJ), p\ is determined by setting fii = [i® in Eq. (|2rS[). 
and ^5 is related to the chiral index as n\ = ((2n + rn)/3) (where (x) represents the 
closest integer to the real number x). Note that we have already ignored the term 
82 § 0^ } I / a(xi,X2) in Eq. (0) by selecting a proper phase for the wave function. The 
coefficient of o\ plays the role of the effective mass of a particle propagating along 
the axis. Therefore, p\ — A® is related to a position-independent mass, and 5m(x2) 
represents a local modulation of the mass. They combine to give the local energy 
gap given in Eq. @. 

It is important to note that the position-dependent mass (or local energy gap) 
corresponds to the existence of a local deformation-induced magnetic field, because if 
a nonvanishing local energy gap remains, then its first derivative in the axis direction 
(^2) is also a nonvanishing value. This corresponds to the magnetic field integrated 
along the x\ direction as 

d 2 Sm(x 2 ) = f^%? b = f j^-B ± ( Xl ,x 2 ). (26) 

The Stokes theorem allows us to imagine a deformation-induced magnetic field Bn 
in a cylinder, defined by the loop integral of the deformation- induced gauge field 
around the axis, given by 

B\\{x 2 )S = j dx l A l {x 1) x 2 ), (27) 

where S denotes the cross-sectional area of the cylinder, assumed to be constant along 
the axis. When the magnetic field escapes from within the cylinder, it penetrates 
the surface. This corresponds to -Bj_ on the surface (penetrating the surface), and 
therefore it changes the mass of the particle (see Fig. |2J) . 

We now briefly consider the interesting case in which p\ — A\ — m{x2) — ^4f m 
changes sign at some position. It is known that an unusual solution, the Jackiw- 
Rebbi modes 14 ) (electron-deformation bound states), exists at the band center. It 
is likely that such a condition will be satisfied for the bulk part of some types of 
peapod carbon nanotubes. 

§4. Examples of Deformation-Induced Gauge Fields 

In this section, we apply the theory to a curved surface of a SWNT. First, we clas- 
sify deformations into two types, using the Slater-Koster scheme. The Slater-Koster 
scheme defines the hopping integral between nearest-neighbor sites as V a (ri) = 
^(flcc + 5r a )e(ri +a ) ■ e(rj), which can include the fact that conducting electrons 
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(a) Energy gap 



+vf \p\ — A\ — 5m{x2)\ 



X2 



-vp \pi — A\ — Sm(x2)\ 



(b) Deformation-induced magnetic field 





Fig. 2. The relationship between (a) the local energy gap and (b) the corresponding deformation- 
induced magnetic fields Bu and B±. It is noted that this is not the usual electro-magnetic 
magnetic field. 



form a 7r-orbital whose wave function extends in the normal direction of the surface. 
Here, V n (r) is a function of the bond length r [K-(a cc ) = V n ] and e(rj) is a unit 
normal vector at fj. We have assumed here that the effect of the cr-bond is included 
in the definition of (r). 




Surface 



Surface 



Fig. 3. A schematic diagram of the deformed surface of a two-dimensional graphite sheet. In these 
diagrams, only two sites (• and o) are depicted. We denote the direct distance between them 
by a cc + 8r a . We refer to (a) as the bond-bending type and (b) as the bond- stretching type in 
this paper. In (a), the direct distance is given by a cc + Eq. 13UII . and in (b) it is not related to 
the normal vector. 

Now, e(ji + R a ) can be expanded in a Taylor series around e(rj) as 



e(n + R a ) = e(r t ) + (R a ■ V)e(r,) + -(R a ■ V)(R a ■ V)e(r t ) + ■■■ . (28) 



Using the normalization condition of the normal vector e(rj) • e(rj) = 1, we obtain 
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e(ri) ■ (R a • V)e(rj) = 0. We can then write 5V a {ri) = V a (ri) — as 

SV a (n) = Sr a (n)(dV v \aJ + V^e( n )(R a ■ V)(R a • V)c(r f ) + • • • , (29) 

where • • • represents higher-order corrections. For nanotubes, the derivative (V) is 
scaled by the inverse of the circumference flC/J) or by the inverse of the length (\T\), 
depending on the direction. We can therefore consider, as an order-of-magnitude 
estimation, the relevant energy of the second term on the right-hand side of Eq. (|29|) 
to be proportional to the inverse square of the circumference. (More specifically, we 
approximate it as the value V n a^ c /\Ch\ 2 ■) This contributes a small (but observable 7 ^ ) 
correction to the energy gap. However, the bond length modulation, denoted by 
3 r a( r i), may give a dominant contribution to the energy gap. 10 ) 

Here, we define two types of deformation depending on the functional form of 
& r a{ r i)- The first type of geometry includes the deformed surfaces of SWNTs in 
which the bond length modulation between adjacent carbon atoms, as shown in 
Fig- Eta), is defined by 



C- &CC 

a ~ 12 



^e(r i )(R a -V)(R a -V)e(r i 



+ •••. (30) 



By substituting this into Eq. (|29|). we obtain 

SV a ( ri ) = D 7T (R a • V)e(ri) • (R a • V)e(n), (31) 
where we have introduced the constant 



Ar = -- 
2 



-^a cc dV n \ acc + Vn 



(32) 



for the sake of convenience. The numerical value of can be estimated from the 
experimental data of Ouyang et al. 7 ) as = V n /8, with V w = 2.60 eV. We 

refer to deformed surfaces whose hopping integral is approximated by Eq. (|31|) with 
an appropriate D n as the bond-bending deformation type. The hopping integral of 
the bond-bending type is calculated from the normal vectors on the surface, which 
means that we only need to parameterize the surface and do not need to perform 
microscopic calculations to obtain the hopping integral. Note we assume here that 
the (indirect) distance between two adjacent atoms defined by the line integral of 
the differential line element on a deformed graphite surface is fixed by a cc and is 
independent of its position. 

In another type of deformation, the dominant contribution to 5V a (ri) comes from 
a change in the bond length. For this type, the hopping integral is approximated by 
the first term on the right-hand side of Eq. ()29|) as 

6V a (r i ) = (dV !T \ a J6r a (r i ) + -- - . (33) 

Equation (|3*3"|) may be relevant to the conducting electrons in a peapod-like struc- 
ture, 10 ) in which the elastic strain is believed to modify the bond length. We refer 
to deformed surfaces whose hopping integral is approximated by Eq. (|33|) as the 
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Table I. Deformation types and expected energy scales. The bond-bending type can be applied 
to ordinary single-wall nanotubes, and the bond-stretching type is believed to be effective for 
peapod nanotubes. 



Type 


SV a (n) 


Target 


Energy scale 


Bond-bending 


Eq. ED 


Nanotubes 


Dndcc/ \Ch 2 


Bond-strctching 


Eq. (E3 


Peapod, etc. 


(a cc dVTr\a cc )Sr a /a cc 



bond- stretching deformation type. The bond length (and its hopping integral) can 
be calculated with a first-principle calculation. 

Table [I] lists the deformation types. A typical energy of D n a 2 c /\Ch\ 2 is ex- 
pected for a deformed Hamiltonian of the bond-bending type, and an energy of 
(a cc dV n \ acc )5r a /a cc for the bond-stretching type. 

4.1. Application to several geometries 

In this subsection, we present examples of the deformation-induced gauge field 
with calculated results for the local energy gap. We calculate the deformation- 
induced gauge field for several geometries, nanotube and peapod, each of which is 
classified as a bond-bending or bond-stretching deformation. We then apply the de- 
formation to the local energy gap formula and compare the results with experimental 
data. We also briefly discuss the effects of the on-site energy. 

4.1.1. Narrow energy gap in a metallic zigzag nanotube (bond-bending type) 

For geometries with a bond-bending deformation, we can apply Eq. ()31j) to the 
deformation- induced gauge field. In this case, by substituting Eq. (|TB|) into Eq. (|3Tj) . 
and using Eq. (f2T|) . we obtain 

v F A x = F™e u ■ e u + F™e u ■ e v + F v x v e v ■ e v , (34) 
v F A y = F™e u ■ e u + F™e u ■ e v + F™e v ■ e v , (35) 

where we have introduced the following quantities: 

i£« = F (^p 2 +pq— ^q 2 ) , Fy U = -FVSq (p + ^ , 

F X V = —F(2pn + qn + pm — qm), Fy V = F\/3(qn + pm + qui), 

F% v == F (n 2 + nm- ^m 2 ^j , F%° = -F^/Hm (n + ^rn^j . (36) 

In the above equations, we define F = 47r 2 Z) 7r /SN 2 and the dimensionless coordinates 
(u, v) through the relations ■ V = |Cft|9i = 2i:d u and T ■ V = |T|(9 2 = 2nd v . 
Finally, e u and e v denote d u e and d v e. The most dominant contributions to the 
deformation- induced gauge field are contributed by the F x u and Fy U terms. This is 
because F itself scales as the inverse square of the nanotube surface area, p and q 
scale as |T|, and n and m scale as \Ch\- Therefore, each component scales as 



rpuu . rpuv . rpvv _ "cc . ~cc . -cc • _ (oj\ 

1 1 1 ~ \C h \ 2 • \C h \\T\ " |T| 2 ' l - x ' y - [6I> 
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This scaling indicates that we can ignore the F™ v and terms if \T\ 3> \Ch\- 

Because the chiral and translational indices determine each component listed in 
Eq. 1|36[) . specifying the normal vector e on a specific surface fixes the deformation- 
induced gauge field completely. We can specify the surface of a tube with a constant 
radius by the vector p(u,v) = (\Ch\/2ir cosu,\Ch\/2nsinu,v). The unit normal 
vector of this surface is calculated as e = (cos u, sinu, 0). We then obtain e u ■ e u = 
1, e u ■ e v = 0, e v ■ e v = 0, yielding the deformation-induced gauge fields vfA x = F£ u 
and vpAy = Fy U . To calculate each component of Eq. (|3l)|) explicitly, we use the 
example of zigzag nanotubes, whose chiral index is (re, 0), with metallic index = 
(where re is a multiple of 3). In this case, Eq. (|3l)|) gives F£ u = —7r 2 D n /n 2 and 
Fy U = 0. Applying this result to the formula for the energy gap of a metallic zigzag 
nanotube, we obtain E gap (v) = 2\F^ U \. We now compare this with the experimental 
result of Ouyang et al. 7 ' The experiment shows that the energy gap for metallic 
zigzag SWNTs with re = 9,12,15 can be fitted by = 7r 2 F 7r /4n 2 , where V n = 2.60 
eV, which corresponds to \D n \ = V n /8. Using the explicit form of the parameter D n , 
we can read off the first derivative of V n (r) as 

mr I l- 9V * if D * <0 ' ^ 
accdKkc = S 1CTr „ „ (38) 



-15K if D n > 0. 



This can be used to analyze the other chiral structures and also for geometries of the 
bond-stretching deformation. It is noted that these values differ significantly from 
that estimated with the scheme employed in Refs. 15) and 16), acc5Kr| acc ~ — 3V^. 
Elucidating the reason for this discrepancy 12 ) is beyond the scope of our low energy 
theory. However, we point out that the overall factor of Fgap may be affected by 
the details of the STS experiments. The important point is that the n dependence 
of the energy gap is the same for the experiment (although there exist data for only 
three different values of n) and theory. Another possible reason for this discrepancy 
is that actual nanotubes are actually between bond-bending and bond-stretching 
classes. We can introduce into Eq. (|32|) a phenomenological parameter a to describe 
such a nanotube as = — | [^acc^Kr|a c + ^r] • Then, we can choose a ~ 3 to get 
(a cc 9K-|a cc ) ~ -3V n . 
4.1.2. Peapod (bond-stretching type) 

We consider the peapod geometry as an example of a bond-stretching deforma- 
tion whose hopping integral depends strongly on the bond length modulation. As 
with some (encapsulated metal fullerene) peapod-like structures, 10 ) elastic strain is 
expected to modify the bond length. We do not attempt to calculate the bond length 
for such systems using a first-principle calculation, but instead we wish to estimate 
the bond length modulation (necessary for explaining an observed local energy gap 
modulation within our theoretical framework). 

We consider zigzag nanotubes (n, 0) of constant radius and assume that A± = A x 
depends only on the axis coordinate. In this case, we can simplify the formula for 
the energy gap into the form E gap (y) = 2vp \p x — A x (y)\. From Eqs. 1)21(1 and 
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the deformation-induced gauge field around the axis can be written 

-v F A x {y) = (dV^laJ (Sn(y) ~ \sr 2 {y) - ~<5r 3 (2/)) . (39) 

It is reasonable to suppose that the bonds, pointing in three different directions, 
are related, because the modulation generates a force between adjacent atoms. For 
zigzag nanotubes, we assume the relation 5r2(y) = Sr 3 (y) = ^5ri(y), which gives an 
equilibrium configuration and preserves the rotational symmetry about the axis. In 
this case, the energy gap formula becomes 



£gap(2/) = 2 



,0 !/„ ^ 5r l{y) 



V FP X ~ -( a ccdV n \ a 



2 a. 



(40) 



where the first term on the right-hand side is given for zigzag nanotubes with chiral 
index (n, 0) by 



o 2^ v / 3//2n\ 2n\ 



(41) 



To realize an observed energy gap modulation on the order of 0.4 eV, 10 ^ the maximum 
bond-length modulation should be on the order of Srf 11 ^ ~ a cc /50, where we have 
used the first case in Eq. (|38|) . Although this quantity was fixed by the experimental 
data of Ouyang et al., 7 - 1 it can be estimated using another scheme 15 ) as a cc dV n \ acc ~ 
—3V n , in which case the maximum bond length modulation should be about o cc /20. 
In any case, the small bond length modulation gives a rather strong contribution to 
the local energy gap. Here, we should note that Cho et al. 17 ) accounted for the local 
energy gap observed by Lee et al. 10 ) in terms of the gap between the NT state and 
the empty metallofullerence state in the gap of the nanotube. 

4.1.3. On-site energy (asymmetry between two sublattices) 

We now briefly comment on the on-site interaction. The on-site interaction is 
defined by 

Wsitc = ^2e(ri)a\ai, (42) 



where e{rt) denotes the on-site energy. The constant component of the on-site energy 
determines the origin of the energy and does not affect the dynamics. By contrast, 
its modulation part can impart physical effects on the conducting electrons. The 
modulation can be divided into a symmetric part and an asymmetric part for the 
two sublattices. They appear in the low-energy Hamiltonian as the coefficients of do 
(the identity matrix) and 03, 18 ) respectively, in the form 

H K ^v F a-(p-A) + G(r)a + P(r)a 3 . (43) 

The former, G(r), acts as the local electrical potential and the latter, P(r), represents 
an asymmetric potential between the two sublattices. While G(r) does not change 
the energy gap of the theory, P(r) does. The local energy gap along the axis can be 
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calculated using dimensional reduction as 




+ 



p 



(44) 



We note that the deformation-induced gauge field can increase or decrease the en- 
ergy gap, depending on the sign of A\. However, the asymmetric potential always 
increases the energy gap. 

§5. Mixing of Fermi points 



To this point, we have disregarded the matrix element that causes mixing of the 
two Fermi points K and K', whose effect is on the order of 5V a {2kp). In this section, 
we set the Fermi velocity, vp, to unity. We begin by introducing the Schrodinger 
equation 



ih di U K ' 



n K o 

o {H K > 



(45) 



where we have defined (J~Lk')' as Hk 1 [see Eq. (|22)l] with the replacement (pi,P2) — ► 
{—pi,p 2 )- Thus, each diagonal block is given, respectively, by TCk = o"i(pi — A\ — 
Af n ) + a 2 (p 2 -A 2 - Af 1 ) and (H K )' = <n(pi -A x + Af*) + a 2 (p 2 + A 2 - A™). 
Roughly speaking, p\ measures the difference between the momenta (about the axis) 
for the K and K' states (Fig. QJ. The u\ terms serve as mass terms (or the energy 
gap) when we consider an effective one-dimensional model. In this regard, we note 
the different signs for an external magnetic field Af 1 . Because Af 11 can be tuned by 
a magnetic field along the axis, it can cause an asymmetry between the masses of 
the two modes; that is, it can break the time-reversal symmetry. 



+4r 




pij—. 


.-^Pi 


— kF 


+k F \ 



Fig. 4. Momentum representation used in § We consider zigzag nanotubes, whose Fermi points 
denoted as solid circles, are located on the k x axis. The dashed lines represent the two energy 
bands near the K and K' points. 



The matrix defined above acts on the four component wave function 



(46) 
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where the two components of ipK, an d V'l?) an d of ifa/, and ip^' , correspond 
to the wave functions at the two sublattices, A and B. The ipR component couples 
to the wave functions on the B-sublattice, ip^ and ipg' , ano - the coupling between 
ip\ and ipg' creates the mixing effect. Therefore, the general dynamics (except for 
the topological defects mentioned below) can be described by the Hamiltonian 



( 


(Wk)i2 





26i \ 


(Wk)21 





26 2 








2b* 2 





(^K')'l2 


V 26^ 





(^K')21 


o / 



(47) 



where 6 = (61,62) are complex functions. 

The mixing of Fermi points may have a strong effect on the dynamics, especially 
when there exist pentagons or heptagons in the surface. 20 ^ We attempt to construct 
an effective low-energy model including the effects of a topological defect and of a 
surface deformation. For this purpose, we first note that the effective dynamics of 
Eq. (|47j) can be expressed as a special case of the Schrodinger equation 

a»D^ = 0, (48) 

where = (p^ — A^)tq — A l ^Ti is the covariant momentum and the summation 
variable // is understood to take the values 0, 1 and 2. Here, A 1 ^ is a generalization of 
the previously defined deformation-induced gauge field, and Tj represents the Pauli 
matrices (with tq the identity element) acting on the wave functions t (ip K ,ip K ). We 
have defined = and po = ihdt, where t is the time variable. The Schrodinger 
equation is formally equivalent to the Weyl equation in U(l) Abelian and SU(2) 
non-Abelian deformation-induced gauge fields. 

The Hamiltonian of Eq. (|47|) can be obtained by setting the generalized gauge 
field of Eq. (gHJ as 

43 _ jem a0 _ Aem a0 _ a a3 _ a 
A \ — ^1 1 ^2 — /1 2 ) A l — Al) ^2 ~~ A 2i 

A\ = -3R(6i + 6 2 ), A\ = 9(6! + 6 2 ), A\ = U(h - 6 2 ), A\ = 9(6x - 6 2 ), (49) 

and setting other components to zero. Concerning the relationship between a pen- 
tagon or heptagon and the non-Abelian deformation-induced gauge field, the main 
contribution from a pentagon or a heptagon is the interchange of the wave functions 
of the two Fermi points. This occurs when a low-energy electron moves around topo- 
logical defect, 20 ^ which can be expressed in terms of a nonvanishing time component 
of the matrix element t\ or r 2 part as 



a\A\- l Al)\(^ 



a\Al + iAl) ;V.V> K ' 



(50) 



This is because a pentagon (heptagon) creates a coupling between ip^ and ip^' , 
and between ip^ and ip^' . Hence, the local dynamics of the low-energy conducting 
electron on a deformed (graphite) surface is governed by Eq. and the lattice 

structure fixes the Abelian and non-Abelian deformation- induced gauge field. 
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§6. Discussion 

The main result for the local energy gap given in Eq. © is determined by 
only the configuration of the deformation-induced gauge field. This is due to the 
special energy dispersion relation of the graphite sheet, which allows us to predict 
several important physical consequences, not by solving the Hamiltonian explicitly 
but from the gauge field configuration itself. To understand the configuration of the 
deformation-induced gauge field, we have classified the geometries into two types: 
a bond-bending deformation and a bond-stretching deformation. For the bond- 
bending type, we can extract the gauge field from the geometric information provided 
by the normal vector as we did in § 14.1.11 where no microscopic theory is necessary 
to calculate the field configuration. By contrast, for the bond-stretching type, we 
need a microscopic theory capable of predicting the bond length or the gauge field 
in order to calculate the local energy gap (§ I4.1.2|) . We can carry out a qualitative 
analysis of the bond-stretching type using a microscopic model and compare it with 
our local energy gap formula. Such a study will be reported elsewhere. 

When we consider the mixing effect due to a short-range interaction, we should 
generalize the geometry-induced gauge field from an Abelian to a non-Abelian field. 
Constructing a general low-energy theory is of prime importance, because a general 
deformation may generate a topological defect, such as a pentagon or heptagon, 
which can mix the wave functions at the two Fermi points. In addition, we can 
consider a higher genus material whose kinematics are nontrivial. 21 )' 22 ) One possible 
way to examine such a material is to extract useful information from a dynamical 
model. 

Finally, we comment on the wave functions in nanotubes with a local energy gap 
modulation. First, note that we did not solve the effective Hamiltonian of Eq. I)22|) . 
but, instead, we extracted the energy gap along the axis by considering the spectrum 
of the effective one-dimensional theory, given by Eq. ((7J), obtained by dimensional 
reduction. Solving the equation for a general A without dimensional reduction would 
be difficult, even in the absence of mixing, because the eigenfunctions of the Hamil- 
tonian, WkV' = Eipi should also satisfy 23 ) the equation 



(p-A) 2 a + hB ± a 3 V = (51) 



where we have set vf = 1 and A cm = 0. The spectrum of the first term on the right- 
hand side is known to be quite nontrivial, even for simple vortex configurations. 24 ) 
Furthermore, it is already difficult to solve the one-dimensional Hamiltonian, due 
to the local modulation mass term. The nature of the wave function should be 
important for investigating the possibility of "energy gap engineering" , 10 ) which 
requires multiple quantum dots prepared by a local modulation of the energy gap. 

§7. Conclusion 



We have examined the effects of surface deformation on the ground state of 
conducting electrons using the nearest-neighbor tight-binding Hamiltonian. Within 
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the framework of the effective-mass theory, we clarified the relationship between a 
local deformation of the lattice and the local energy gap along the axis in terms of the 
deformation- induced magnetic field. We formulated an effective theory describing 
the dynamics on a general deformed surface, including topological defects that can 
cause a mixing of the two Fermi points. The theory is formally equivalent to the 
Weyl equation in C/(l) Abelian and SU{2) non-Abelian deformation-induced gauge 
fields. 
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